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ABSTRACT 

We estimated the dynamical surface mass density (S) at the solar Galactocentric 
distance between 2 and 4 kpc from the Galactic plane, as inferred from the observed 
kinematics of the thick disk. We find S(2;=2 kpc)=57.6ib5.8 Mq pc~^, and it shows 
only a tiny increase in the z-range considered by our investigation. We compared our 
results with the expectations for the visible mass, adopting the most recent estimates in 
the literature for contributions of the Galactic stellar disk and interstellar medium, and 
proposed models of the dark matter distribution. Our results match the expectation 
for the visible mass alone, never differing from it by more than 0.8 Mq pc~^ at any 
z, and thus we find little evidence for any dark component. We assume that the dark 
halo could be undetectable with our method, but the dark disk, recently proposed 
as a natural expectation of the ACDM models, should be detected. Given the good 
agreement with the visible mass alone, models including a dark disk are less likely, but 
within errors its existence cannot be excluded. In any case, these results put constraints 
on its properties: thinner models (scale height lower than 4 kpc) reconcile better with 
our results and, for any scale height, the lower-density models are preferred. We believe 
that successfully predicting the stellar thick disk properties and a dark disk in agreement 
with our observations could be a challenging theoretical task. 



Dipartimento di Astronomia, Universita di Padova, Vicolo Osservatorio 3, 1-35122, Padova, Italia 



- 2 - 



Subject headings: dark matter — Galaxy: kinematics and dynamics — Galaxy: struc- 
ture 



1. INTRODUCTION 



Today it is widely accepted that dark matter is a fundamental component of the Universe, 
which plays a key role in the processes of galaxy formation and evolution. Cosmological A'^-body 
simulati ons accurately pre dict the evolution of the dark component and its actual spatial distri- 
bution (jMoore et al.lll999l ). In the A cold dark matter (ACDM) cosmology, spiral galaxies accrete 
smaller building blocks into both their spheroidal compone nt and their stellar disk. The presence 
of an old thick disk, very common among spiral galaxies (IDalcanton fc BernsteinI 120021'). is often 
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considered the product of one or more merging events (e.g. 
proposed that, as satellites are torn apart by tidal forces, the y should deposit t heir dark matter 
into a flat dark structure. This idea was recently explored bv iRead et al.l (j2008l . hereafter Re08), 
who showed that the presence of a dark disk is a natural expectation of the ACDM model. As 
a result of their simulations, they proposed that a galaxy such as the Milky Way should host a 
relatively thin dark feature (exponential scale height 2.1-2.4 kpc), with a local density at the solar 
position (pd) 0.25-1.0 times that of the dark halo (ph)- The proposed dark component is much more 
flattened than the dark halo, but it is still noticeably thicker than any visible disk, because the 
scale height of the Galactic old Thick stellar disk is ~0.9 kpc, while young er stars and inter stellar 



medium (ISM) form even thinner structures (0.3 and 0.1 kpc, respectively, lJuric et al.ll2009l ). 



In the last year it became progressively accepted that, if the ACDM cosi nology is the c orrect 
model, dark disks should be ubiquitous in spiral galaxies. More recently, iPurcell et al.l (120091 , 
hereafter Pu09), elaborated new models, proposing a thicker (scale height 4.6 kpc) but less dense 
dark disk in the Milky Way {pd/ph <0.30). 

The presence of the dark disk is strongly related to the formation of the stellar thick disk. 
However, its merging origin is currently u nder debate, because o f the difficulties of the models in 
reproducing all its properties. For example. iBournaud et al.l (|2009l ) argued that thick disk formation 



throu gh turbulent and clumpy phases at high redshift can explain its lack of flaring (jMomany et al 



2006 ). at variance with the merger scenario. This model would not necessarily require the presence 
of a dar k disk. On the other hand, the presence of a phantom disk is also an expectation of MOND 
theory (iMilgromI Il983l ) , where the departure of gravitation from Newtonian law should cause the 
detection of an additional amount of disk matter (jMilgroml l200ll ) . It is thus clear that the Milky 
Way dark disk has become a benchmark for many theories, from gravitational law to cosmological 
galaxy formation, and thick disk origin. 

We are performing ai i extensive survey to study the kinematical and chemical ve rtical structure 
of the Galactic thick disk ( Carraro et al hood ). Preliminary results were presented by lMoni Bidin et al, 
(|2009l ) . In this Letter we analyze the vertical trend of the surface mass density as inferred by thick 
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disk kinematics, in search of evidence for any dark component. The detailed analysis of the kine- 
matical results will be published in a later paper (C. Moni Bidin et al. 2011a, in preparation, 
hereafter Paper I), and the full study of the hypothesis and equations presented here will follow 
(C. Moni Bidin et al. 2011b, in preparation, hereafter Paper II). 



OBSERVATIONAL DATA 



Our investigation is based on a sample of ~1200 red giants defined by iGirard et al.l (|2006l ). 
vertically distributed with respect to the Galactic plane in a cone of 15 degrees radius centered 
on the south Galactic pole. All objects have Two Micron All Sky Survey (2 MASS) photometr y 
(jSkrutskie et al.ll2006l ) and absolute proper motions from the SPM3 catalog (jGirard et al.l 12004 ). 
The sample was selected by applying a color cut in the infrared color-magnitude diagram, to iso- 
late intermediate-metallicity thick disk stars. Main-sequence (MS) dwarfs were excluded both by 
a sloped cut at fainter magnitudes, which excludes all but the nearest {d <63 pc) dwarfs, and by 
conservative kin ematical criteria im posing a stellar velocity lower than the local escape velocity 
(550 km s^^; see lGirard et al.ll2006l . for more details). We collected high-resolution Echelle spectra 



for two-thirds of the Girard et al.'s sample, during 38 nights at La Silla and Las Campanas obser- 
vatories. The distribution of proper motions and colors of this sub-sample was analyzed, to ensure 
that no selection effect was introduced. We visually inspected all the spectra and excluded the 
residual low-metallicity stars ([Fe/II]< —1.5), most probably halo contaminators, and misclassified 
dwarfs. Radial velocities (RVs) were measured for all the stars with a typical error of 0.5-0.7 km s~^, 
cross -correlating each spectrum with standard stars observed during the same runs (jTonrv &: Davis 
1979 ). Distances were estimated with a color-absolute magnitude relation calibrated on the red 
giant branch of 47 Tuc, whose stellar po pulation is very similar to the intermediate-metallicity 
Galactic thick disk (jWyse &: Gilmordl2005i ). Finally, the three velocity components in the Galactic 
cylindrical coordinate system {U,V,W), and their associated errors, were calculated for each star 
from its proper moti on, RV, and dista nce, and the uncertainties on these quantities. The error in 
distance was ~ 20% ( Moni Bidinll2009l ). while the err or in proper motio ns was fixed to 3 mas yr~^ 
(T. M. Girard, 2009, private communication, see also lGirard et al.ll2006l . for a detailed discussion). 
The mean errors in (U,V,W) thus increased from ~(7,36,38) km s~^ at z=2 kpc, to ~(10,61,66) km 
s~ ^ at 4 kpc. The d etails of observations, data reduction, and RV measurements were presented 
bv lMoni BidinI (j2009l ). and they will be fully discussed in Paper I. 



In this Letter, we restrict our analysis to the ~300 stars beyond 2 kpc from the Galactic plane, 
to avoid contamination by the thin disk. We considered only stars with \W\ <150, \U\ <300, and 
—500 < V <300 km s~^, to exclude halo stars and/or bad measurements. We did not apply a cut 
in velocity errors, because it systematically excludes high-velocity stars, biasing the results: the 
propagation of the uncertainty in distance in fact introduces a term proportional to the velocities 
themselves. The sample was binned with respect to distance from the Galactic plane z, following 
three different criteria: three large bins of 85 stars each, five bins of 50, and five of 45 stars. The 
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Fig. 1. — Dispersions of the three velocity components, as a function of Galactic height, and their 
least-squares fit. Full dots and line: ajj] open squares and dashed line: cry; open stars and dotted 
line: aw- The errors in z (<100 pc at 4 kpc), given by the statistical error on the mean, are omitted 
for clarity. 
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dispersion in the three velocity compon ents was calculated in each bin , fitting the correspon ding 



probability plot with a linear relation (|Lutz &: HansonI Il992l : see also iBochanski et alj 120071 . for 



an application to a very similar astrophysical case). In the nearest bins {z <2.3 kpc) the fit was 
performed outside the ztlcr range, to avoid the residual thin disk contamination (5%-10% according 
to our estimates). At increasing distance an overestimate was a serious possibility because of some 
outliers, most probably stars with an incorrect distance and/or proper motion. We consequently 
excluded from the fit the points on the wings of the distribution showing clear deviation from 
linearity, i.e. outliers departing from the underlying normal distribution. The mean error in each bin 
was quadratically subtracted to derive the final estimate of the dispersion. Finally, the dispersion 
in each velocity component was plotted against the average z of the bin, and a least-squares linear 
relation was fitted to derive the vertical trend of dispersions. 

The results were very similar in the three binning criteria: differences on the derived quantities 
were smaller than 1.5o", and in most cases agreed within la. We therefore gathered the information 
of all the bins together, fitting a final plot comprising 13 points. The results are shown in Figure [H 
We derive a small vertical gradient of ajj, ay, and aw of 6.2ib3.7, 4.5±1.9, and 2.8ib0.5 km s~^ 
kpc~^, respectively. The statistical uncertainty of the least square solution was assumed as the 
error on the linear profile parameters. These could be underestimated because of the correlation 
of the fitted points, but they are only about 30% higher when fitting the uncorrelated bins of 50 
stars, and this has only minimal impact on the final results, because the uncertainties on the thick 
disk parameters dominate the error budget. 

Around z=2.5 kpc we find a sudden deviation from linearity of the cross-term of the dispersion 
matrix, UW. This feature could be due to a stellar sub-structure, such as a comoving group of 
stars. The use of only 13 bins is inadequate to reveal the general trend, hence we binned the data 
in overlapping groups of 50 stars at steps of 2 stars, as shown in Figure [2j We finally derived the 
linear relation required for our calculations excluding the "depression" at 2.5 kpc. The line is a 
good fit outside this feature. 



3. THE THEORY 

Our formulation is based on the following assumptions: 

I Steady state. The thick disk is in equilibrium with the Galactic potential, as expected for an 
old stellar population. Therefore, all temporal derivatives are set to zero. 

II Locally flat rotation curve. The rotation curve is assumed flat at the solar Galactocentric 
distance. 



Ill No net radial or vertical stellar flux. The mean radial an d vertical velocity cq i nponents are 
zero, while the rotational component shows a non-null lag (|Chiba &: Beersll2000l : iGirard et al. 
2004 ). 
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IV Exponential radial dispersion profiles. All the velocity dispersions decrease with R following 
an exponential law, with a scale length h/j^o- equal to the one observed for the mass density 
(hi?,p). 

V Vertical constancy of scale lengths, h^ p and h^ do not depend on z at the sola r position. 
This is observationally verified for the mass density (ICabrera-Lavers et alJl2005l ). and it is 
assumed valid for the velocity dispersions, because their radial behavior is linked to the mass 
distribution by assumption pvp . 



VI Null cross term on the Galactic plane. UW{z = 0)=0. This hypothesis is required for 



symmetry reasons (see e.g. lBienaymell2009l ). and it is observationally confirmed by lFuchs et al 

kood ). 



The hypothesis (jlVI) is observationally confirmed for aw (|van der Kruit &: Searld Il982l ). Its ex- 
tension to the other components relies on the controversial assumption of constant anisotropy, i.e. 
= 0. We have no information abo ut the radial trend of t he velocity dispersions, but this 



assumptio n is supported by observatio ns (iLewis &: Freemanlll989l ). and numerical simulations for 
R <9 kpc dCuddeford fc AmendtlflQQ^ V 



Integrating the Poisson equation from —z to z, assuming that the vertical component of the 
force is null on the plane, we obtain: 



27rGE(z) 



km^^^''^'^'^^'^'^ 



(1) 



where is the surface mass density between itz, and F/j is the radial component of the force. In- 
serting the Jeans equations in Equation ([1]), making use of our hypotheses (jIt)- (|Vip . and calculating 
simple derivatives and integrals, we eventually obtain the final expression: 
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where hz p is the thick disk exponential scale height, and: 
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To calculate the surface density we must insert into Equation ^ the vertical trends of the kine- 
matical quantities au, cry, cr\Y, and UW, plus the three parameters Rq, h^_p, and h^^p. We fixed 
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R0=8.O=bO.3 kpc, and we defined the th ick disk scale height and length as the mean of about 20 lit- 
erature estimates (see lMoni Bidinll2009l . and Paper II for the bibliographical references), obtaining 
h_R,p=3.8±0.2 kpc and h2^p=0.90±0.08 kpc. The quoted errors are given by the error on the mean, 
that should be the most appropriate statistic when averaging many uncorrelated measurements. 
However, they could be too small because of the traditional uncertainty on these parameters, al- 
though the estimates converged considerably in the last years. We therefore also considered an 
uncertainty of 0.4 and 0.12 kpc, respectively, comprising all the measurements of the last decade 
within ±2(7 (~70% within itlcr). These larger errors necessarily decrease the significance of the 
results, but without altering the general conclusions, as can be deduced from Figure [H The error 
on T,(z) was calculated from the propagation of errors of all the quantities in Equation 

In deriving Equation ([2]), we assumed the kinematical quantities as symmetric with respect to 
the plane, so that their integrals between —z and z are twice the product of integration between 
and z. This is easily justified for the dispersions, but it may not be the case for UW. We therefore 
also considered an antisymmetric cross-term, where the third term of Equation ([2]) vanishes. As 
shown in Figure [3l a symmetric cross-term fails to return physically meaningful results, because it 
violates two minimum requirements: the surface density must at least account for the known visible 
matter, and it cannot decrease with z. We will therefore assume the cross-term as antisymmetric, 
and consider only the results obtained under this hypothesis hereafter. 



4. RESULTS 



Our results are shown in Figure H We find S(2 kpc)=57.6±5.8 Mq pc^^ and the curve IS 
nearly fiat in the whole range, while the error increases constantly to ~12 Mq pc~^ at 4 kpc. The 
solution increases 1 Mq pc~^ between 2 and 3.2 kpc, then it turns slightly downwards. However, 
the decrease is so small (0.55 Mq pc^^) that a small refinement of any parameter would correct it, 
for example increasing h^^p by 0.1 kpc. Moreover, this problem is present only in the lower half of 
the family of solutions defined by the ilo" strip. We did not alter the input parameters to amend 
it, because this would add a high degree of arbitrariness to the results. 

In Figure [3] we overplot a model of surface mass density due to the visible mass. The thin and 
thick disk expon ential scale heights , the halo shape, and the local thick disk and halo normalization 
were t aken fromlJuric et al.l (120091 ). We included a thin layer (100 pc) of ISM contributing 13 Mq 



pc 



-2 



(|Holmberg fc Flvnnl boOof T" and th e thin disk density on th e plane was normalized by the 



requirement Ildisk(l-1 kpc)=40 Mq pc~^ (jHolmberg fc Flvnrj|2004l ) . This quantity includes all disk 
stellar compo nents and remnaiits, and it is currently the best estimate often assumed in Galactic 
mass models (jPehnen fcBinneJl998l : loiling k MerrifieldlboOll : Iweber de Boerlbood ). 



The agreement between the visible mass and our dynamical solution is striking, and there is 
no need to invoke any da rk component. In Figure [3] we also plot a model including a dark halo 



(jOlling &: Merrifieldll200ll ). The curve is completely incompatible with our results, both for its high 



- 8 - 



value and steep slope. One can argue that we are actually deducing the mass from the change in 
the gravitation potential: the dark halo could be too uniform and extended to cause a detectable 
change on stellar kinematics in a range of 4 kpc. This point requires a detailed analysis that is 
beyond the scope of this Letter, and we leave it to a later study (Paper II). Here we assume that 
the dark halo is undetectable and focus on the dark disk, for which this explanation is not viable. 

Figure H] compares our results with models where a dark disk is considered. They require 
the definition of two additional parameters: the exponential scale height of the dark component 
hz^D, and its density at z=0, normalized with respect to the dark halo, Pd/ph- We considered two 
families of models, a relatively thinner (h2^£)=2.4 kpc) dark disk with prf/p/j= 0^25- 1.0, as proposed 
by Re8, and a thicker, less dense one (h2^£)=4.6 kpc, p^/ Ph=0 -15-0. 30), from |Pu09l . The measured 
surface density matches the baryonic mass alone, and any curve including the dark disk departs 
from the central solution, being therefore less likely. Hence, there is no evidence for a flattened dark 
component, but its expected contribution is small compared to the errors and its existence cannot 
be completely ruled gut. H owever, not all combinations of dark disk parameters are per mitted by 
our observations: the |Pu09| solutions fall in the Icr-strip, but only the less dense models of lReOSi do. 
This is entir ely due to the lower densities of the first family of models. In fact, for a fixed Pd/ph, the 
thick |Pu09l model s depa rt more from the observed curve, and they require a higher increment of 
E(z) compared to iReOSi models (see Figure [5]). We can thus derive important constraints: a lower 
density is favored in all cases, and thinner disks should be preferred, while thick and dense dark disks 
(Mz,D ^4 kpc, Pd/Ph ^0.5) are less likely. Numerical simulations showed that low-latitude merging 
of m a,ssive satellites are require d to form a heated disk kinematically similar to the Galactic thick 
disk (IVillalobos &: Helmill2009l ). but these events produce denser dark disks (|Pu09l ). In addition, 
ReOSi showed that the observed density of the stellar thick disk can be reproduced only by a few 
models, all implying a more massive dark disk {pd/ph >0.4). In summary, correctly predicting 
the thick disk kinematical and dark disk propertie s is a challenging task . Models for thick disk 
formation alternative to the merging sce i iario (e.g. jBournaud et al.ll2009l ) are preferred, while the 
MONDian prediction of lBienavme et al.l (|2009l ). i.e. a 60% increase of disk mass due to a phantom 
disk (~30 Mq pc-2 at 2 kpc) is in contradiction with our results. 

It could be argued that our model of visible mass, it relies on poorly-constrained quantities, 
such as the ISM contribution and Sdisk(l-1 kpc), and a downward correction of these parameters 
would shift all the model curves to lower values. Therefore, in Figure [5] we analyzed both the 
expected and measured increment of S(z), i.e. the surface density of the mass comprised between 
z and 2 kpc. This quantity is completely independent of the ISM at these Galactic heights, and 
the assumed value of 5](iisk(l-l kpc) introduces a negligible uncertainty (<0.15 Mq pc~^). Figure[5] 
shows exactly the same situation analyzed previously, and our conclusions should be considered 
independent of the details and uncertainties of the modeled visible mass. 
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5. CONCLUSIONS 

We estimated the Galactic surface mass density between 2 and 4 kpc from the plane, finding 
E(2 kpc)=58±6 Mq pc~^, and a nearly flat curve. Our results strikingly match the visible mass 
alone, and we do not detect evidence for any dark component, although the dark halo could have 
passed unseen. There is no compelling evidence for a dark disk, but within the errors of our 
investigation, its existence cannot be completely excluded. We derive important constraints on its 
expected properties: lower densities {pd/Ph <0.25) should be preferred in any case, and a thin dark 
disk {hz^D <2.5 kpc) better reconciles with observations. A thick and dense dark disk (hz^o >4 kpc, 
Pd/Ph >0.5) should be excluded. Any merging model aiming to reproduce the formation of the 
Galactic thick disk and a fiat dark component will need to consider the constraints from our 
investigation. 
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Fig. 2. — Measured cross-term UW as a function of distance from Galactic plane. The line indicates 
the least squares fit obtained omitting the depression at 2.5-3 kpc. 
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Fig. 3. — The surface density vertical profile derived from our observations. The full line shows 
the results under the assumption of cross-term antisymmetry, while the dash-dotted line indicates 
the results for a symmetric cross-term. The dashed and dotted lines indicate the corresponding la 
strips. Light gray curves reprensents the expectations of the visible mass alone (lower curve), and 
of a visible mass+dark halo model (upper curve) . 
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Fig. 4. — Comparison of our results with the expectations of dark disk models. The black full line 
shows the observational result with its lo" errors (dashed lines) as in Figure El The dotted lines 
indicate the la strip when the enhanced errors on the thick disk parameters, discussed 
in the text, are considered. The light grey curv es are models where a dark disk is added to 
the visible mass, basing on the models propo sed by iReOg (dashed lines, hz^£)=2.4 kpc and three 
different values of local density), and bv |Pu09l (full line, h.z,D=^-^ kpc, two values of local density). 
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Fig. 5. — Same as Figure HI but for the increment of the surface mass density as a function of 
with respect to 2 kpc, instead of its absolute value. 



